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Abstract 

We consider improving the sensitivity of future interferometric gravitational-wave detectors by 
simultaneously injecting two squeezed vacuums (light), filtered through a resonant Fabry- Perot 
cavity, into the dark port of the interferometer. The same scheme with single squeezed vacuum 
was first proposed and analyzed by Corbitt et al. Here we show that the extra squeezed 

vacuum, together with an additional homodyne detection suggested previously by one of the au- 



thors 



allows reduction of quantum noise over the entire detection band. To motivate future 
implementations, we take into account a realistic technical noise budget for Advanced LIGO (Ad- 
vLIGO) and numerically optimize the parameters of both the filter and the interferometer for 
detecting gravitational-wave signals from two important astrophysics sources, namely Neutron- 
Star-Neutron-Star (NSNS) binaries and Bursts. Assuming the optical loss of the ^ 30m filter 
cavity to be lOppm per bounce and lOdB squeezing injection, the corresponding quantum noise 
with optimal parameters lowers by a factor of 10 at high frequencies and goes below the technical 
noise at low and intermediate frequencies. 
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I. INTRODUCTION 



During the last decade, several laser interferometric gravitational- wave detectors includ- 
ing LIGO la], VIRGO 4], GEO600 Q and TAMA ^\ have been built and operated almost at 
their design sensitivity, aiming at extracting gravitational-wave signals from various astro- 
jhysical sources. At present, developments of next-generation detectors such as AdvLIGO 
7j are also under way, and sensitivities of these advanced detectors are anticipated to be 
limited by quantum noise nearly over the whole observational band from 10 Hz to 10^ Hz. 
At high frequencies, the dominant quantum noise is photon shot noise, caused by phase 
fluctuation of the optical field; while at low frequencies, the radiation-pressure noise, due 
to the amplitude fluctuation, dominates and it exerts a noisy random force on the probe 
masses. These two noises, if uncorrelated, will impose a lower bound on the noise spectrum, 
which is called the Standard Quantum Limit (SQL). In terms of gravitational- wave strain 
h = AL/L, it is given by 



S' 



SQL 



8h 



(1) 



It can also be derived from the fact that position measurements of the free test mass do not 
commute with themselves at different time Q]. 

The existence of SQL was first realized by Braginsky in 1960's jo], [l^. Since then, vari- 
ous approaches are proposed to beat the SQL. One recognized by Braginsky is to measure 
conserved quantities of the probe masses (also called quantum nondemolition quantities . 



This can be achieved, e.g. by adopting speed-meter configurations 
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16|, 



which measure the conserved quantity - momentum rather than the position. An alterna- 



tive is to change the dynamics of the probe mass, e.g. using optical rigidity 17|, 



m 



which case t 
Chen 
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le free mass SQL mentioned is no longer relevant. As shown by Buonanno and 
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211 ]. optical rigidity exists in signal-recycled (SR) interferometric gravitational- 



wave detectors; therefore we can beat the SQL without radical redesigns of existing topol- 
ogy of the interferometers. Another approach is to modify input and/or output optics 
of the interferometers such that photon shot noise and radiation-pressure noise are corre- 
ated. After the initial paper by Unruh 221], this was further developed by other authors 
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28|, |29|, |30|] . A natural way to achieve this is injecting squeezed vacuum. 



whose phase and amplitude fluctuations are correlated, into the dark port of the interferom 



eters. With great advancements in preparation of the squeezed state 
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FIG. 1: Schematic plot of the proposed configuration. Two squeezed vacuums s,p are injected 
from both side of the filter cavity rather than the one which was considered in Refs. . The 

signal is detected by the main homodyne detector (MHD) and additional detector (AHD) is made 
in the idle port of the filter cavity. 

interferometers will be promising candidates for thi rd-g eneration gravitational-wave detec- 



tors. As elaborated in the work of Kimble et al. [27|, frequency- dependent squeezing is 
essential to reduce the quantum noise at various frequencies of the observation band. In 
addition, they demonstrated that this can be realized by filtering the frequency-independent 
squeezed vacuum through two detuned Fabry-Perot cavities before sending into the inter- 
ferometer. Their results were extended by Buonanno and Chen 
general cases were discussed. 
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where the filters for 



Another method, which also uses an additional filter cavity and squeezed light, but in 
a completely different way, was proposed by Corbitt, Mavalvala, and Whitcomb (hereafter 
referred as CMW) They proposed to use a tuned optical cavity as a high-pass filter 
for the squeezed vacuum. This scheme does not create the noises correlation, but instead, 
renders the noises spectral densities frequency dependent. At high frequencies, the phase 
squeezed vacuum gets reflected by the filter and enters the interferometer such that high- 
frequency shot noise is reduced; while at low frequencies, ordinary vacuum transmits through 
the filter and enters the interferometer, thus low-frequency radiation-pressure noise remains 
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unchanged. One significant advantage is that the squeezed vacuum does not really enter the 
filter cavity and thus it is less susceptible to the optical losses. However, it does not perform 
so well as hoped for, and there is a noticeable degradation of sensitivity in the intermediate- 
frequency range. One of us - Khalili pointed out that this has to do with the quantum 
entanglement between the optical fields at two ports of the filter cavity. Equivalently, it 
can be interpreted physically as the following: some information about phase and amplitude 
fluctuation flows out from the idle port of the filter cavity and the remaining quantum state 
which enters the interferometer is not pure. In order to recover the sensitivity, the filter 
cavity needs to have a low optical loss such that these information can be collected with an 
additional homodyne detector (AHD) at the idle port. Given an achievable optical loss of the 
filter cavity ~ lOppm per bounce, Khalili showed that we can obtain the desired sensitivity at 
intermediate frequencies. A natural extension of this scheme is sending additional squeezed 
vacuum into the idle port of the filter cavity such that the low-frequencies radiation-pressure 
noise is also suppressed. The corresponding configuration is shown schematically in Fig. [1], 
where two squeezed vacuums s and p are injected from two ports of the filter cavity, and some 
ordinary vacuum n leaks into the filter due to optical losses. By optimizing the squeezing 
angles of squeezed states, we will show that the resulting quantum noise is reduced over the 
entire observational band. 

The outline of this paper is as follows. In Sec. [Tll we will calculate the quantum noises 
in this double squeezed-input CMW scheme with AHD (later referred as CMWA). We will 
use the same notation as in Ref 21|, which enables us to extend the results in Ref. ^ 
to the case of signal-recycled interferometers easily. In Sec. IIIIl we numerically optimize 
the parameters of this new scheme for searching the gravitational-wave signals from NSNS 
binaries and Bursts. Finally, we will summarize our results in Sec. IIVI For simplicity, we 
will neglect the optical losses inside the main interferometer, but we consider the losses from 
the filter cavity and also non-unity quantum efficiency of the photodiodes. The losses from 



the main interferometer are not expected to be important as shown in Refs. 19|, Um, l33 |. 
The main notations used in this paper are listed in Table [B 
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TABLE I: Main notations used in this paper. 



Quantity 



Value for Estimates Descriptions 



n 

c 

(Jo 

m 
L 

Ic 

ic = 

Tarm 

<5 

7 



8(^0 Je 

mLc 



L 



f 



3.0 X 10^ m/s 

1.8 X 10^^ s-i 

40 kg 

4 km 

840 kW 

(27r X 100)3s~3 

27r X 100 s-i 



30 m 



7l,E,L - 
7/ = 71 + 7E + 7L 

n {i = s,p) (lnl0)/2 (lOdB) 

9i [i = s,p) 



Gravitational wave (sideband) frequency 

Speed of light 

Optical pumping frequency 

Mass of the end mirrors 

Length of the arm cavity 

Circulating optical power 

Bandwidth of the arm cavity 
Amplitude reflectivity of the SRM 
Phase detuning of the SR cavity 

Effective frequency detuning of the SR interferometer 
Effective bandwidth of the SR interferometer 
Homodyne angle of the MHD 
Homodyne angle of the AHD 
Length of the filter cavity 

Bandwidth of the filter cavity 

Squeezing factors 
Squeezing angles 



II. QUANTUM NOISE CALCULATION 



A. Filter cavity 

In this section, we will derive single-sided spectral densities of the two outgoing fields d 
and q as shown in Fig. [H From the continuity of optical fields, we can relate them to ingoing 
fields, which include two squeezed vacuums s, p and one ordinary vacuum h entering from 
the lossy mirror (LM) Specifically, we have 

a{n) = ni{n)s{n) + r{n)p{n) + A{n)h{n), (2) 
q{n) = nE{n)p{n) + T{n)s{n) + AE{^)h{n). (3) 

Here d = (a^, O'lp)'^, Q = {qa, (l<fi)'^ are amplitude and phase quadratures (subscript A 
stands for amplitude and (p for phase). In our case, the carrier light is resonant inside the 
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filter cavity. Therefore, the effective amplitude reflectivity TZ, transmissivity T and loss A 
can be written as, 

-7? m^ 7i - 7e - 7l + ^ 7e - 7i - 7l + '^^ ^ 
^i(^) = 7^ , 7^E(^^) = , (4a) 

nn) = (4b) 
A,(n)^Z^, Mn)-'-^ (4c) 

where 7/ = 71 + 7e + 71- They satisfy the following identities, 

|7^,(^])|2 + |T(fi)p + = |7^E(fi)|' + |T(fi)|2 + |^E(^^)r = l, (5a) 
n;{Q)T{Q) + 7^E(^^)T*(^]) + ^*(1])^e(^^) = o. (sb) 

If the input and end mirrors of the filter cavity are identical, namely 71 = 7e, we will 
have TZ\,TZe ~ 1 and T ~ when 1] ^ 7/ and TZ\,TZe ~ and T ~ 1 when Q <^ 7/. 
Therefore, the squeezed vacuum s enters the interferometer at high frequencies while p 
becomes significant mostly at low frequencies. By adjusting the squeezing factor and angle 
of these two squeezed fields, we can reduce both the high-frequency shot noise and low- 
frequency radiation-pressure noise simultaneously. 

To calculate the noise spectral densities, we assume these two squeezed vacuums have 
frequency-independent squeezing angles 6i {i = s, p) and can be represented as follows: 

s = R{rs,9s)vs, p = R{rp,ep)vp, (6) 

with 

~ / cosh r — cos 9 sinh r — sin 9 sinh r \ 

R{r, 9)^{ . (7) 

\ — sin 9 sinh r cosh r + cos 9 sinh r / 

Here r is the squeezing factor, and Vi = {viA, Vi^)^ are ordinary vacuums with single-sided 
spectral densities: S'A(f^) = S^{Q) = 1, SAipi^) = j^]. 

The noise spectral densities of the two filter-cavity outputs can then be written as 

Saai^) = |7^,(^])|2i?(2r„2e,) + \Tin)\^R{2rp,29p) + |^(f])pT, (8) 
Sqq{n) = \nEm'Ri2rp,29p) + \T{n)\'R{2r,,29,) + lAem'l (9) 

Sag{n) = 7^E(^^)T*(^])^(2rp, 29p) + n;{n)T{n)R{2rs, 29^) + AE{yt)A;{n)i, (lo) 
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where Saq{^) is the cross correlation between two outputs; I is the identity matrix and 

^ ( ^^"<"' ] (11) 

with i = aa, qq, aq, whose elements are single-sided spectral densities. 



B. Quantum noise of the interferometer 



According to Ref. 2l|], the input-output relation, which connects ingoing fields a and 
gravitational- wave signal h with the outgoing fields b, for a signal- recycled interferometer 
can be written as 

b = lf(ca + B-!^]. (12) 



In the above equation, 

M = [5^ - {n + i-ff]n^ - 5 Lc, (13) 

and C is the transfer function matrix with elements 

Cn = C22 = - 52 + 7') + 5 ic. (14a) 

C'l2 = -27 5fi^ = 275^2 -27 i„ (14b) 

where ic = SuQlc/imcL). The elements of the transfer function vector D are 

Di = -26y/Yi~M, D2 = -2(7 - in)y^M- (15a) 

The effective detuning 5 and bandwidth 7 are given by 

^ _ 2 rsR 7arm sin(2 (Psr) . 
~ l + riR + 2rsRCos(20sR)' ^ 

l + ^2^ + 2rsRCOs(2 0sR)' ^ ^ 

where rsR is the amplitude refiectivity of the signal recycling mirror (SRM) and 0sr is phase 
detuning of the SR cavity. If the main homodyne detector (MHD) measures 



b^{n) = 07 [sin bA{ft) +cos(f) b^{Q)] + ^/l-riv{Q), (17) 

where is the homodyne angle and v is the additional vacuum due to non-unity quantum 
efficiency of the photodiode, and then the corresponding /i-referred quantum-noise spectral 



density can be written as 

^ (sin0 cos (f))C Si sin (f) cos^)^ + i^lMP 

Shin) = /i^QL— — -i-^ (18) 

(sin0 cos (j)) D D^sin (j) cos0)'^ 

It can be minimized by adjusting the squeezing angle 6 of s and p. We can estimate the 
optimal 6 qualitatively from the asymptotic behavior of the resulting noise spectrum. At 
very high frequencies (fi 7), from Eq. ([2]), d ~ s and thus 

Sh{n) oc cosh(2r,) + cos[2(0 + 9^)] sinh(2r,). (19) 

If the squeezing angle of s 

= I + nvr - (20) 

where n is integer, we achieve the optimal case, namely Sh oc e~^^^ Similarly, at very low 
frequencies (fi ^ 7), we have d ~ p and the spectral density Sh oc e"^^*" if 

2 cos sin P 



9p = arctan 



(21) 



cos{(3 — 0) + 3 cos(/3 + 

More accurate values for optimal Og^p can be obtained numerically as we will apply in the next 
section. Given optimal 6s,p, the sensitivity of this double squeezed-input scheme improves 
at both high and low frequencies. However, due to the same reason as in the case of single 
squeezed-input that two outputs of the filter cavity are entangled this double squeezed- 
input scheme does not perform well in the intermediate frequency range. To recover the 
sensitivity, we need to use an additional homodyne detector (AHD) at the idle port E of the 
filter cavity. The corresponding measured quantity is 



q^ifl) = [sin C bA{n) + cos ( b^{n)] + ^/T^ i'{n) (22) 

where ( is homodyne angle and v'{Q) is the addition vacuum, which enters due to non- unity 
quantum efficiency of photodiode. We combine g^(f2) with the output 6</,(f2) using a linear 
filter )C{Q), obtaining 

d{n) = b^{n) - ic{n) q^{n). (23) 

Corresponding, the noise spectral density of this new output d{Q) can be written as 

Sain) = Sf,^iQ) - 23?[/C(n)^s.,a.(^)] + \m)\'S^,in). (24) 



8 



The minimum quantum noise is achieved when )C{Q) = Sf^^ / Si^{Q) and the resulting 
/i-referred noise spectrum with AHD is then 

^AHD^^^ _ ^^^^^ _ ^2 r/|(sin0 cos^^y g|E(sinC cosQ^P ^25) 

SQL^^.^^ COS 0) D-D''" (sin cos^)"^ 'S'c(^) 

where S(;{fl) = rj{sm( cosC)5'E(sinC cosC)""" + 1 — 77. The second term has a minus sign, 
which shows exphcitly that the sensitivity increases as a result of additional detection. 



III. NUMERICAL OPTIMIZATIONS 

In this section, we will take into account realistic technical noise and numerically optimize 
interferometer parameters for detecting gravitational- wave signals from specific astrophysics 
sources which include Neutron-Star-Neutron-Star (NSNS) binaries and Bursts. 



For a binaries system, according to Ref. [3J], spectral density of gravitational- wave signal 
is given by 

^^(2-/) = i^UW^- ^^^^ 

Here the "chirp" mass A4 is defined as 7V1 = /i^^^M^/^ with /i and M being the reduced mass 
and total mass of the binaries system. With other parameters being fixed, the corresponding 
spectrum shows a frequency dependence of f~'^^^ ■ Therefore, as a measure of the detector 
sensitivity, we can defined an integrated signal-to-noise ratio (SNR) for NSNS binaries as 

The upper limit of the integral /max ~ /isco ~ 4400 x {M/Mq) Hz is determined by the 
innermost stable circular orbit (ISCO) frequency, and the lower limit /min is set to be 10 Hz, 
at which the noise can no longer be considered as stationary. Here we choose M = 2.8Mq, 



the same as in Ref. 351]. Here S')^"'^'^* is the quantum noise spectrum derived in the previous 



sections and Sf^^ corresponds to the technical noise obtained from Bench 



3- 



Another interesting astrophysical sources are Bursts [36| . The exact spectrum is not well 
modeled and a usual applied simple model is to assume a logarithmic-fiat signal spectrum, 
i.e. Sh{2iif) oc f~^. The corresponding integrated SNR is then given by 

PBurst. 0^ / ^quant^27r/) + 5*-h(27r/)' ^ ' 



mill 
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The integration limit is taken to be the same as the NSNS binaries case. 

To estimate the SNR and also motivate future implementation of this scheme, we assume 
the filter cavity has a length of ~ 30 m and an achievable optical loss lOppm per bounce and 
also consider the non-unity quantum efficiency of the photodiodes t] = 0.9 for both MHD 
and AHD. Other relevant parameters will be further optimized numerically. For compari- 
son, we will also optimize other related configurations, which includs AdvLlGO, AdvLlGO 
with frequency-independent squeezed-input (FISAdvLIGO for short), and CMW scheme. 
Specifically, free parameters for these different schemes that need to be optimized are the 
following: 

AdvLlGO: rsR, (/-sr , , (29a) 

FISAdvLIGO: rsR,0sR,0,^, (29b) 

CMW: rsR, 4>SR, (p, lu IE, 0,,ep, (29c) 

CMWA: rsR, 0sr , , 7i, 7e, , , C • (29d) 



TABLE II: Optimization results for different schemes 



Configurat ions 


NSNS 


Bursts 


P 




4>SR 


71 


7E 


</> 


9p 


c 


P 


rsR 


(t>SR 


71 


7E 


<^ 


Op 


c 


AdvLlGO 


1.0 


0.8 


1.4 






-1.0 






1.0 


0.7 


1.5 






-0.2 






FISAdvLIGO 


1.0 


0.7 


1.5 






-1.2 


-0.6§ 




1.5 


0.8 


1.6 






0.0 


-1.6§ 




CMW 


1.0 


0.8 


L6 


240 


240 


-0.8 


-0.6 




1.5 


0.8 


1.6 


0.0 


0.0 


0.0 


-1.4 




CMWA 


1.2 


0.7 


L6 


230 


210 


0.0 


-0.1 


0.7 


1.5 


0.8 


1.6 


140 


140 


0.0 


-0.2 


0.9 


^ This is the 


squeezin 


g an 


gle e 


in the 


case of 


sing 


e squeezed-input. 



The resulting optimal parameters for different schemes are listed in Table [ITl They 
are rounded to have two significant digits at most to balance various uncertainties in the 
technical noise. The integrated SNR p is normalized with respect to that of the AdvLlGO 
configuration. The optimal 6s from the numerical result is in accord with the asymptotic 
estimation, namely 6's ~ (vr/2) — (c.f. Eq. fl20|) ). 

The corresponding quantum-noise spectrums of different schemes optimized for detecting 
gravitational waves from NSNS binaries are shown in Fig. [2J The FISAdvLIGO, FISAd- 
vLIGO and CMW schemes almost have the identical integrated sensitivity and CMWA 
scheme shows moderate 20% improvement in SNR. This is attributable to the fact that 
the signal spectrum of NSNS binaries has a /"^^^ dependence and low-frequency sensitivity 
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FIG. 2: Quantum-noise spectrums of different schemes with optimized parameters for detecting 
gravitational waves from NSNS binaries. The optimal values for the parameters are listed in Table. 

IE 
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FIG. 3: Quantum-noise spectrums of different schemes which are optimized for detecting gravita- 
tional waves from Bursts. The optimal values for the parameters are listed in Table. [Ill 

is very crucial. However, due to low-frequency technical noise, advantages of the CMWA 
scheme are at most limited in the case for detecting low-frequency sources. 

The case for detecting gravitational waves from Bursts is shown in Fig. [21 AH other 
three schemes have a significant 50% improvement in terms of SNR over AdvLIGO. The 
sensitivities of the optimal FISAdvLIGO, CMW and CMWA at high frequencies almost 
overlap each other. In addition, detuned phase 0sr of the signal recycling cavity of those 
three are nearly equal to 7r/2, which significantly increases the effective detection bandwidth 
of the gravitational- wave detectors and is the same as the resonant-sideband extraction 
scheme. This is because a broadband sensitivity is preferable in the case of Bursts which 
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FIG. 4: Quantum-noise spectrums of different schemes using the same parameters as the optimal 
CMWA to show how different parameters affect sensitivity of the CMWA scheme. 

have a logarithmic flat spectrum. 

To show explicitly how different parameters affect sensitivity of the CMWA scheme, 
we present quantum-noise spectrums of different schemes using the same parameters as 
the optimal CMWA in Fig. |H In the case of AdvLIGO, we obtain a resonant-sideband 
extraction (RSE) configuration with 0sr ~ 7r/2. The quantum noise of FIS AdvLIGO with 
squeezing angle 6 = 6s is lower at high frequencies but higher at low frequencies than the 
RSE AdvLIGO. FISAdvLIGO with 9 = Op behaves in the opposite way with significant 
increase of sensitivity at low frequencies but worse at high frequencies. The CMW scheme 
with double squeezed- input, just as expected, can improve the sensitivity at both high and 
low frequencies but performs not so well at the intermediate frequencies. The CMWA scheme 
performs very nice over the whole observational band compared with others and it would 
be more attractive if the technical noise of the AdvLIGO design can be further decreased. 



IV. CONCLUSION 



We have proposed and analyzed the double squeezed-input CMWA scheme as an option 
for increasing sensitivity of future advanced gravitational- wave detectors. Given an achiev- 
able optical loss of the filter cavity and lOdB squeezing, this CMWA configuration shows 
a noticeable reduction in quantum noise at both high and low frequencies compared with 
other schemes. Since the length of the filter cavity considered here is around 30 m, with the 
developments of low-loss coating and better squeezing-state sources, it could be a promising 



12 



and relatively simple add-on to AdvLIGO without a need for dramatically modifying the 
existing topology. 
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